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abstract 

We study diffractive effects in two dimensional polygonal billiards. We derive an analytical trace 
formula accounting for the role of non-classical diffractive orbits in the quantum spectrum. As an 
illustration the method is applied to a triangular billiard. 
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During the last decade several methods based on periodic orbit (PO) theory have been successfully 
employed to study quantum systems whose classical equivalent is chaotic (see e.g. ||]). PO theory 
applies also when the system is not fully hyperbolic (when some orbits appear in families Q) or 
integrable ||. More recently it has been refined to include complex orbits Q and diffractive effects 
[||, |6| . In this line we aim at studying the problem of wedge diffraction as an extension of the standard 
PO theory. This is one of the oldest and simplest example of diffraction (see e.g. [Q]) and it is also 
the case where the diffractive corrections to semiclassics are the more important. 

In this Letter we calculate for the first time the role of non-classical diffractive orbits in the spec- 
trum of two-dimensional polygonal billiards. We derive a trace formula embodying the contribution 
of diffractive PO's to the level density (Eq. (0)). This contribution is of order Vh smaller than the 
contribution of isolated PO's and is the next order term in the trace formula. As an example the for- 
malism is applied to a triangular billiard with angles (vr/4, ir/6, 7ir /12) and one sees that it provides 
a very accurate description of the Fourier transform of the spectrum. 

We consider a quantum particle enclosed in a polygonal billiard B and we impose Dirichlet boun- 
dary conditions on the frontier dB. Hence the associated Green function is solution of the following 
equation : 



(A B + /c 2 )G(q s ,q A , k) = <5(q s - q A ) inside B, (1) 
G(q s ,q A ,/c) = on dB . 

where q is a coordinate in configuration space. 

The semiclassical approximation for G reads (see e.g. ||) 

e i(kL — htt/2) 

G {q B ,q A ,k)= / . , — ' ( 2 ) 

qA ^ qs ivSmkL 

where the sum is taken over all classical trajectories going from q A to q fl . In (||) L is the length of 
the trajectory and fj, is the associated Maslov index P]. In polygonal enclosures the boundary has 
no focussing components, there are no caustics and /x is simply twice the number of bounces of the 
trajectory on dB. 

In polygonal billiards the hamiltonian flow is discontinuous on the vertices ||] and when the angle 
at a vertex is not of the form it/n (n £ IN*) this causes diffraction (see e.g. [p!o|). Then, following 
Keller's geometrical theory of diffraction fll|| , one is lead to consider non-classical contributions to 
the Green function which are "diffractive orbits" starting at q A , going to a vertex qi and then to 
q B . These orbits are non classical because at qi the reflection is not specular. Far from the region of 
discontinuity of the hamiltonian flow the corresponding Green function may be taken to be : 



G!(c[ B ,q A ,k) =G (qi,q A , k)V x (0, 0') G o(qB, qi, k) , (3) 

where T>i(9,9') is a diffraction coefficient evaluated in the solvable case of two semi infinite straight 
lines meeting with an angle 7 equal to the interior angle of the polygon at qi. 9 (resp. 9') is the angle 
of the incoming (resp. outcoming) trajectory at q\ with the boundary. T>i(9,9') reads [ll], 12| : 
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4 S in(vr/iV) S m(g/jV) S in(g7jV) 

^(MH-^- — * eT¥~. — n o=F-' (4) 

(cos _ _ cos — )(«» - - COS — ) 

where N = 7/V is not assumed to be an integer. 

As stated above one sees in expression (||) that when 7 is of the form ir/n, T>\ is zero and there 
is no diffraction. Indeed in this case a trajectory passing by qi is the limit of a trajectory bouncing 
specularly n times near the vertex and a contribution of type (Q) accounts for the effect of the wedge. 
This is to be related to the fact that in this case there exists a n th iterate of the flow which is continuous 
||. Note also that T>\ is zero if 6 or 0' is equal to or 7 (i.e. in the case of a diffractive trajectory 
having a segment lying on a face). 

For an orbit with several diffractive reflections at points qi...q^ formula (||) becomes 

G u {n B ,q A ,k) = G (qi,cu,A;) j Yi VjGofa+iiVjityj A/G (qB,cu,,&) > ( 5 ) 

where T>j is the diffraction coefficient at point q., as given by (Q). 

In (|2|), (|3|) and (||) the indices 0, 1 or v of the Green function recall that diffractive effects are 
subdominant (by a factor of order k~ u l 2 ). There might be less severe non analyticities on the boundary 
leading to higher order diffractive corrections. Note also that we are using here a simple approximation 
for the Green function which is not valid when the angles 6 and 9' at an edge are such that the diffractive 
orbit is close to be real ; in this case the coefficient T>i(6,9') diverges. This occurs in vicinity of the 
line of discontinuity of the hamiltonian flow. In order to have a formula valid in all regions of space 



one should use a uniform approximation such as first provided by Pauli [12] and whose general form 
is given in jl0| (see also |l3]| ). 

The level density p{k) is then obtained from the Green function by the usual formula : 

Ik r 

p(k) = Im d 2 gG(q,q,fc). (6) 

vr JB 

p(k) can be separated in a smooth function of k, p(k) plus an oscillating part p(k). The zero 
length trajectories in (^) contribute to p and will not be considered in detail here (see [0]). When 
G is replaced by its semiclassical approximation (|2|) a stationary phase evaluation of (^) corresponds 
in considering only the contribution of classical PO's to p. When diffractive orbits such as (||) are 
taken into account one is lead to consider also "diffractive PO's" which are PO's with one or several 
diffractive reflections (example of such orbits are given on Fig. 1). 

Let us consider first the contribution of classical PO's. In a polygonal enclosure there is a drastic 
difference between PO's with even and odd number of bounces. The latter ones do not remain periodic 
when a point of reflection is translated along a face (they period-double into a PO with twice as many 
bounces). This can be understood by remembering that, for the phase-space coordinates transverse 
to the direction of an orbit, a bounce on a straight segment leads to an inversion. On the other hand, 
PO's with an even number of bounces form families which correspond to local translation parallel to 
the faces of the polygon. They are neutral (or direct parabolic, see [Q) PO's to which the usual trace 
formula does not apply ; we use a generalization of Gutzwiller theory which is valid for the case of 
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degenerate PO's 0. We quote here the result and leave detailed discussion for the future [13]. A 
family of orbits contributes to p(k) as : 



/ kL 

pik) <- J ^3 d± cos(krL - vr/4) . (7) 

(0) is written for the general case of the r th iterate of a primitive orbit of length L(r£ IN*). d± is 
the length occupied by the family perpendicular to the orbit's direction. It is equal to dcoscj), where 
d is the length occupied by the family on a face and 4> the angle between the direction of the orbit 
and the normal to this face. 

For an isolated PO with an odd number of bounces one has the following contribution : 

p( k ) < ^cos{krL). (8) 

Formula @ holds when the number of repetitions is odd. When r is even, the r th iterate of an 
isolated orbit leads to a family and formula ([?]) applies. 

The derivation of the contribution of a diffractive PO is patterned on what is done in Gutzwiller's 
trace formula for an isolated PO. The length of a closed diffractive orbit in the vicinity of the diffractive 
PO is expanded up to second order and the trace of the Green function is evaluated by a stationary 
phase approximation. The final contribution of a generic diffractive PO with v diffractive reflections 
to the oscillating part of the level density reads : 

cos(/cL - pir/2 - 3i/tt/4) . (9) 

In © L\...L v are the lengths along the orbit between two diffractive reflections. L\ + ... + L v — L 
is the total length of the diffractive PO. p is the Maslov index which is here twice the number of 
specular reflections. Formula @ is the most important result of this paper. Note that different 
diffractive orbits may combine if they have diffraction points in common. Hence repetitions of a 
primitive diffractive orbit appear as a special case of @ ; in this case however, in the first factor L/ir 
of the r.h.s. of (||), L should be understood as the primitive length of the orbit. The above formulae 
show that the contribution of a family of orbits in of order 0{k 1 / 2 ), for an isolated orbit it is 0(1) 
and for a diffractive PO it is 0(k~ v l 2 ). Nevertheless we will see in the following that diffractive orbits 
have a very noticeable contribution to the level density. 

We will now illustrate our approach by studying a specific example. Let us consider a triangle with 
angles (7r/4, 7r/6, 7-7r/12). As explained above diffraction occurs only at the vertex with angle 77T/T2. 
The scale of lengths and wave- vectors is fixed by the value h of the height going from this vertex to 
the opposite face. We take h = 1 in the following. The shortest classical and diffractive PO's in this 
triangle are shown on Fig. 1. Diffractive reflections are indicated with a black spot. Note that the first 
orbits are diffractive, classical orbits (isolated or in families) occur at greater lengths. The spectrum 
was computed numerically by expanding the wave function around the vertices with angles n /A and 
7r/6 in "partial waves" which are Bessel functions with a sinusoidal dependance on the angle defined 
near the vertex considered. More precisely if r n and ip n are polar coordinates defined near the vertex 
ir/n (n = 4 or 6) the partial waves in this region are of the form J nm {kr n ) sm{nnup n ) with m G IN*. 
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One then imposes matching of the wave function and of its first derivative along the height h (see 
details in |L3|]). We determined the first 957 levels, up to k max ~ 96. The accuracy of the computation 
was tested by varying the number of matching points and of partial waves. We evaluate the typical 
error on an eigenvalue as being of order of a hundredth of the mean level spacing. 

In order to visualize the importance of classical and diffractive PO's of successive lengths in the 
spectrum we study the regularized Fourier transform of the level density : 

F(L) = / k e lkL ~ ak p{k)dk . (10) 
Jo 

If kmax — ► +oo and if the regularizing coefficient a is set to zero in ( |l0| ) F(L) is just a series of 
delta peaks centered on the lengths of the classical and diffractive PO's. The multiplicative factor k 
in ( JlOj ) is meant to cancel the singularity k~ u / 2 of the contribution of a diffractive PO of type (|9|) with 
up to v = 2 diffractive reflections. We take here a = 9/k^^ and plot on Fig. 2. The numerical 

result is represented by a thin line and the semiclassical approach (0,§) corrected by diffractive PO's 
(||) by a thick line. We also included the contribution of p(k) in order to reproduce the initial peak 
at L = 0. We see that the agreement is excellent. Note that the existence of diffractive PO's is of 
great importance for reproducing all the peaks in |_F(L)|. This is illustrated on the figure where their 
contribution (||) has been dashed. 

Here several comments are in order. Note first that the diffractive PO's labelled 2 and 4 on Fig. 
1 have not been included because their diffraction coefficient is zero. Also the orbit labelled 7 on Fig. 
1 has a non standard contribution : it is an isolated orbit which accounts for boundary effects on the 
family with the same length (labelled 6 on Fig. 1). In addition to the orbits of this family it has 
an extra reflection on the bottom face (the same type of orbit was considered in Ref. [15, 16]). The 
weight of PO number 7 is reduced by a factor 1/2 compared to @ since one integrates only over closed 
orbits on one side of this limiting PO. Also we included repetitions of diffractive PO's number 1 and 
3 and they can be seen to have still a noticeable contribution. We did not include the diffractive PO 
composed by the sum of orbit 1 and 3 although it can be considered as a small diffractive correction to 
the contribution of family 6. Indeed the orbit "1+3" lies just on the region separating real orbits from 
diffractive ones and as mentioned above it can not be accounted for by a simple diffraction coefficient 
such as (||). This type of corrections will be treated in a forthcoming publication [13]. 



To summarize let us emphasize the important role of non classical orbits in the spectrum of 
quantum billiards. The existence of these orbits affects qualitatively the Fourier transform of the 
spectrum. The above example is only one among others were the discontinuity of the classical dynamics 
is linked to strong diffractive corrections to semiclassics. It was argued in [17| that the same type 



of corrections should be taken into account for the 3 dimensional icosahedral billiard. We expect 
also diffractive effects - of the same order as those described here - in more general billiards with 
cusps (non-polygonal or with an additional external field) ; in these cases a simple generalization of 
formula (0) accounts for the role of diffractive PO's. We note finally that the present work illustrates 
that semiclassical methods provide a very appealing tool which, when corrected with tunneling or 
diffractive effects, allows to describe accurately the solution of partial differential equations using 
simple geometrical methods. 

It is a pleasure to thank E. Bogomolny and D. Ullmo for fruitful discussions. Division de Physique 
Theorique is a Unite de Recherche des Universites de Paris XI et Paris VI associee au CNRS. 
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Figure captions 



Figure 1. The shortest classical and diffractive PO's in the triangle (vr/4, ir/6, 77r/12). All these 
orbits are self-retracing. For diffractive PO's the diffraction point is marked with a black spot. Orbits 
6 and 10 form families, 5 and 7 are isolated. The lengths are given in unit of the height of the triangle. 

Figure 2. 1-^(^)1 as a function of L. The thin line is the numerical result and the thick line 
the semiclassical approximation (|7|,||) with diffractive corrections (|9|). The two curves are hardly 
distinguishable. The contribution of the diffractive PO's has been dashed. 
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